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current {C)ao in the jet from primary parton Aq which takes into account intermediate partonic pro- 
cesses in the sub-jet production where Cp < {C)ao < at high energies. We generalize previous 
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results. 
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1 Introduction 



The coUimation of hadrons in jets is a basic phenomenon of high energy colhsions and its quantitative 
understanding is an important task for QCD. Simple differential characteristics of a jet are their energy 
and multiplicity angular profiles. The coUimation of energy and multiplicity in the jet follows from the 
dominance of gluon bremsstrahlung processes in the parton cascade evolution. Whereas the former is 
more sensitive to the hard processes inside a jet the second one is sensitive to the soft parton emissions 
from the primary parton. 

The characteristics of soft particle production, such as particle multiplicities, inclusive distributions and 
correlation functions, are derived in QCD in the Modified Leading Logarithmic Approximation (MLLA) 
(for review, see UJ) which takes into account the leading double logarithmic terms and the single logarith- 
mic corrections. These azimuthally averaged quantities can be obtained from an evolution equation for 
the generating functional of the parton cascade. This equation provides also Next-to-MLLA corrections 
taking into account energy conservation of parton splittings with increased accuracy. The corresponding 
hadronic observables can be obtained using the concept of Local Parton Hadron Duality (LPHD) IJl 
which has turned out a successful description of many hadronic phenomena (see, e.g. [3]). 

A characteristic prediction of QCD is the increased mean particle multiplicity in a gluon jet over the 
quark jet, asymptotically by the ratio of color factors r = Nq/Nq Nc/Cp = 9/4 modified by 
corrections in powers of ||4j |5] 0I- These results have been compared with LEP data on 2- and 
3-jet events, interpreted as qq and qqg primary production [7,8, 9|. In applications with primary hadron 
beams, such as at HERA, RHIC, TEVATRON or LHC it is necessary to select jets by an angular cut to 
remove the effects from the hadron remnants. Such an analysis of multiplicities has been performed at 
the TEVATRON 1 10] with different jet opening angles and energies E. The multiplicity as well as 
differential multiplicity distributions in MLLA depend on the maximum transverse momentum in the jet 
("virtuality"), i.e. on the combination 

Q^EQ (1) 

in the small angle approximation. This scaling behavior has been verified ifTOll for the multiplicities 
Ng, Nq in the angular range Q = 0.28 . . . 0.47 and also for the peak position of the inclusive rapidity 
spectra [11|. Results from LEP and TEVATRON support in general the theoretical treatment based on 
the MLLA evolution equation and the agreement with data improves with the inclusion of higher order 
terms in 

For small opening angles G the definition of the jet axis becomes problematic because of the fluctuations 
at low particle multiplicity. The multiplicity at smaller angles Q can be defined in a more stable way 
without reference to a jet axis by the double inclusive correlation density between two particles. The 
particles (p2) within the cone G around particle (pi) are counted with angle G12 (including particle (pi) 
itself), then the superposition of the particles (pi) in the jet weighted by their energies Ei is performed. 
For a jet of full energy E, in a form differential in the energy fraction X2, the multiplicity in the cone G 
is given as (see Ref. O) 



dx2 E Jq J J aidEidE2d@i2 \ E J 

. ^ /-.^. 

In this definition, the collection of particles (pi) serves to replace the jet axis and therefore, also small 
angles are well defined, which is not the case for the usual global jet axis definition. 

The evaluation of Eq. ^ for the distribution in energy X2 and, after integration, for the multiplicity of 
particles (p2) in the cone has been obtained in LLA for the integral over the energy Ei in Refs. fT2l [T]|. 
The definition of jet axis is especially relevant for the inclusive distribution of the transverse momentum 
Ut of particles and this quantity has been derived with reference to the energy flux in the MLLA |[T3]| and 
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including an important class of Next-to-MLLA (NMLLA) contributions as obtained from the evolution 
equation 1 14|. The predictions on the distributions have been found in a rather good agreement with 
the preliminary data {Ut > 1 GeV) from the CDF collaboration ifTSll using so far a conventional jet 
definition. In this paper we study the effects from the definition Eq. Q and derive results for particles 
confined into small cones. This allows the study of particle distributions with smaller in the jet and 
lower effective energies (virtuality) E@ than the nominal jet energy E. 

For decreasing cone size © of a sub-jet in the definition of Eq. ^ the simple scaling behavior as in 
Eq. ^ will be broken. The occurrence of intermediate processes leads to a mixing of quark and gluon 
sources and, asymptotically in the LLA, the hadron multiplicity iVj^^ in the sub-jet with opening angle 
Q in the primary jet from parton Aq {Q or G) and large opening angle ©o ~ 1 is given by [ 1] 

Nl^{@,E,eo) = ■^JC)a,-N^{E@), (3) 

proportional to the multiplicity Nq of an isolated gluon jet of virtuality £^0. The average color current 
{C)ao inside the small cone is obtained in the same approximation fl^ as 

{C)Ao = {u)%-N, + {u)'i^.CF, (4) 

where (w)^^ denotes the mean energy fraction of parton A inside parton jet Aq for given jet virtualities 
EQ and EQq. This quantity is found in between the color factors 

Cf < {C)ao < N,. (5) 

The hadron multiplicity N^, as suggested by LPHD, is taken proportional to the parton multiplicity N 
at cut-off scale Qq, i.e. N^{EQ) = K ■ N{EQ, Qq) with hadronization constant K. 

Using the alternative approach for jet definition from Eq. Q we derive in this paper results on the particle 
multiplicity as a simple observable showing its dependence on the color current and its variation with jet 
energy and opening angle © in the improved accuracy of MLLA and NMLLA. Some phenomenological 
applications will be studied. 



2 Evolution equations for multiplicities in quark and gluon jets 

In MLLA the evolution of azimuthally averaged quantities with jet energy E and jet opening © is given 
by an evolution equation for the generating functional for the parton momenta in the jet llj. The evolution 
involves q^, the running coupling constant of QCD 

as = as(Ee) = ^^TWTT> /3o = TTT f —^c - -Tr] (6) 

where A = Aqcd is the intrinsic scale of QCD and /3o is the first term in the perturbative expansion 
of the /3— function H, Nc is the number of colors, T/j = nj/2, where nj is the number of light quark 
flavors. In Double Logarithmic Approximation (DLA) ag is also linked with the anomalous dimension 
7o of twist-2 operators by 

-,i^-i(Ee) = 2N,:idm = ^_ Ye=ln^, A = l„^. (7) 

where Qq is the coUinear cut-off parameter for kx = EQ > Qq. The results depend on energy and 
angle only through the variable Yq, i.e. the maximum transverse momentum in the jet. We also set 
Yq = Yq + Xin the following. 



' Following common practice, we consider in this work on soft particles the 1-loop expression for the running coupling 
constant, therefore A has to be considered a phenomenological parameter not related to the usual A-jyg-. 
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For the evolution of the mean multiplicities in quark and gluon jets one obtains in MLLA the coupled 
system of two evolution equations UJ 



d 



dYe 



NhiYe) = l^'dxji{Y^)[^^{x)[N^{Ye + \nx)+N^{Ye + \na-x))-N^{Ye) 
+ n^$g(x) (N^iYe + Inx) + N^iYe + ln(l - x)) - N^iYe) 



(8) 



^N^iYe) = jyxjiiY^)[<^^ix)[N^iYe + lnx)+N^iYe + Hl-x))-N^iYe: 



Together with initial conditions for the multiplicities at threshold (every jet contains only one parton, 
quark or gluon) these equations determine the multiplicities at the higher values of EQ. The DGLAP 
splitting functions denoted by are given by 



$g(x) = i-(l-x)[2-x(l-x)], ^Q(x) = ^^[x' + {l-x)% $§ = ^(1-1 + 1 

The equations above for multiplicities can be solved analytically 111 in terms of an expansion in 70 of 
the ratio of multiplicities r and the QCD anomalous dimension 7 defined by 

Nh dlnNh 

The Eqs. ^ are complete only up to MLLA order. The emerging higher order terms in the asymptotic 
expansion, although not complete, provide, however, an important contribution improving the constraint 
from energy conservation. In particular, they build the correct threshold behavior of the parton cascade 
and therefore a better behavior at the present energies is obtained. A full numerical solution of the 
evolution equation |[T6l corresponding to an all order resummation has in fact provided a rather close 
reproduction of the multiplicity ratio r above; the same is also true for the result of the HERWIG MC 
ifTTl at parton level and at hadron level. On the other hand, the truncated series in 70 while asymptotically 
convergent, diverges at low energies. 

We limit ourselves here to the NMLLA expression of r that reads 

r = ro(l-ri7o-r27o') + 0(7o) (10) 

where the asymptotic value of r is tq = Nc/Cp = 9/4. The MLLA term ri has been calculated in H. 
The coefficients can be obtained from the Taylor expansions of N^(YQ+lnx) and Nl^(YQ+ln{l—x)) 
for Yq ^ Inx and Yq » ln(l — x) respectively and from equating the terms of the same order in 70 
in both sides of Eq. dUl; the values for Nf = 3 for NMLLA read n = 0.185 and r2 = 0.426 The 
NMLLA solution for the mean multiplicity in a gluon jet is found as S 

^oiYe) - K' exp (^^V^" ^p^) ^''^ 

with ai = 0.28, 02 = 0.38 at Uf = 3; K is the LPHD normalization factor. The pre-exponential term 
{Yq)~"'^/^° is the MLLA contribution to Nq, while the one oc 02, the NMLLA one. We give accordingly, 
the 0(70) expressions of the first and second logarithmic derivatives of Nq that follow from Eq. ([TT]) 

= joiEQ) - a^jiiEQ) + 0(70^), (12) 
^' + =7i{Ee) + 0{jl) (13) 



dY% V dYe J dY% 



which we will all use in the following. The NMLLA solution for Nq can be obtained by substituting Eq. 
([IB into Eq. OO]). 
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3 Energy-multiplicity correlations 



We consider the production of a jet in a high energy collision (pp, pp, ep, e~^e . . .) initiated from a 
parton Aq of energy E and with opening angle Go which separates it from other jets. A sub-jet is defined 




Figure 1: Inclusive production of hadron /i in a sub-jet of opening angle G inside a high energy jet of 
total opening angle Go- The "jet axis" of parton A corresponds to the direction of the energy flux, i.e. 
the energy weighted direction of particles in the jet according to Eq. Q. 



by the opening angle Q < Qq with respect to parton A with momentum fraction u. The distribution 
F^^ {x, Q, E, Go) of the momentum fraction x of a hadron h in this sub-jet is obtained by the integration 
of the double-inclusive correlation function in Eq. |2l see Ref. [12|. This process is depicted in Fig. [T] 
Notations and kinematics are identical to the ones used in the work of Ref. |[T3l which also contains 
further details. 

The distribution F^^ is obtained as the convolution of two fragmentation functions HI 

F^^ix,e,E,eo)= Yl / duDi^{u,Eeo,uEe)D'x(-,uEe,Qo) . (14) 

This equation expresses the correlation between the energy flux and one particle (/i) within the sub-jet 
with angle Q generated from intermediate parton A and with energy fraction u. The function in 
([14)) describes the probability to emit parton A of virtuality uEQ with energy fraction u off the initial 
parton Aq (Q or G), taking into account the evolution of the jet between Q = EQq and Q = EQ. The 
function D\ in Eq. (fT4l) describes the probability to produce the hadron h off A with energy fraction 
x/u and transverse momentum scale k± « uEQ > Qq. 

Integrating (fT4l) over the energy fraction x yields the corresponding sub-jet multiplicity of hadrons 
inside the angular range G < Go of the jet Aq 

iV^„(G;i?,Go) ~ V C duuDiju,EQo,uEe)N'l{uEe,Qo), (15) 

where A^^ is the number of hadrons (partons at scale Qq) produced inside the sub-jet A of total virtuality 

uEQ > Qq. 

We also note two limits of the sub-jet multiplicity in Eq. ([T5] ). At first, in the limit of large angle G ^ Go 
one finds, see also iJl, 

G^Go: dX^6{1-u)6X, N^^^iQ; E,Qq) ^ nX,{EQq;Qq) (16) 

and the sub-jet multiplicity coincides with the overall jet multiplicity. At small opening angle G — > ^ 

G^O: D'X^dil-xMd'X, iVX(G;£;,Go)^l. (17) 

the sub-jet multiplicity approaches the limit of one particle. This limit cannot be reached in our approxi- 
mation with a truncated expansion in 70. 
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3.1 Approximations with the leading parton 



The convolution integral in Eq. (fTSl ) is dominated by the region n 1 as N^{uEQ, Qo) is exponentially 

u~^p{lnu) leads to a slowly varying function uD^ (u) (for 



rising with u and the behavior D^^ (u) ~ 
further discussion, see L13J). Therefore, we use the DGLAP expressions for {u, EQo,uEQ) in the 
following, then the evolution between the two scales depends only on the variable 



1 



4iVc/3o 



In 



In 



1 



In 



uEe 



In 



leo + A 
In ti + Ye + A 



(18) 



and for L>^^ (n, EQo,uE@) we also write D^^ [u, £,{u)) for short. 

Within the leading parton approximation we expand the multiplicity N'^{uEQ, Qq) in ([TS] ) at u ~ 1. As 
in the MLLA evolution equation ^ and thereafter we take a logarithmic dependence of multiplicity A'^^ 
and for Inn ^ Ye = ln(S0/(5o) and for £'0 S> A, this quantity can be written as 



N^{lnu + Ye) ~ N^{Ye)+lnu 



dYe 



dYi 



+ 0(a 



3/2 N 



(19) 



up to the NMLLA level of accuracy. Therefore, in this approximation the correlation integral in ([T5]) can 
be replaced by a sequence of factorized terms showing alternating sign 



Here depends only on Ye and the mean energy fraction (^i)^^ is given by 

where D^^ {j, ^(1)) is the representation of DGLAP fragmentation functions in Mellin's space HI 



Dlij = 2,e(l)) = C duu^-^DX{u,i{l)) 
Jo 



3=2 



The other moments can be obtained by differentiation after the exponent j 



i=2 



(20) 



(21) 



(22) 



(23) 



i=2 



The lower bound in Eqs. (l22l ) and (1231 ) has been set to "0" in the limit Qq <^ EQ, as the threshold region 
in the integral Eq. ([T5] ) is exponentially suppressed. The second 0{as) term in Eq. (|2TI) follows from 
the scaling violation in DGLAP fragmentation functions when one sets n = 1 in the third argument of 
D^^ in Eq. ([TS] ): it has been proved in App. D of fW ] not to exceed 5% of the leading contribution, that 
is why we neglect it hereafter and simply set (u)^^ = {j = 2, ^(1)) as in ifTlfT?!. 
Now Eq. (l20l ) can be conveniently rewritten in the form 



dYe ' 2 

and, using Eqs. (fT2l) and ([T3l) . we obtain an expansion in 70 ~ 

A 

1 



1 d^N'l 
^°iVj^Y|" 



NliYe) 



(24) 



(25) 
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where the contribution oc {u)^^ is the leading one (LLA) determined in m 

(tx)g = /3(l-e-^«), (u)g = a + /3e-^«, (26) 
{u)^ = a(l-e-^«), {u)^ = P + ae-"'^ (27) 

with 7 = |Cf + |n/, a = (3 = where a + /? = 1 and ^ = ^(1). The functions (nlnu)^^^ 
and (uln^ n)^^ are determined from (1231 ). as in f\3l . Taking Nc and Cj? at their QCD values the results 
still depend on ny. As this dependence is rather weak, we present the results with numerical coefficients 
only for n/ = 3 in the appendix. In the limit — > Gq, — > 0) the expression in Eq. (|25T ) fulfills the 
limit Eq. ([T6l ) as the mean values in Eq. (1231 ) vanish and (n)^ — > 6"^ . 



3.2 Multiplicity and color current 

The expression for can be written in compact form as in Eq. ^ introducing the color current. In 
the sum over A in Eq. (1251 ) when A is a quark (^4 = Q) we replace 



Q 



1 + ri7o + r27i)N^, h = rf + ra = 0.46, 



(28) 



according to Eq. (ITOl ) and after multiplication with the terms inside the square bracket we keep only 
terms up to 0(79) so as to respect the NMLLA scheme (1 + 70 + 7q). One thus obtains the average 
multiplicity N^^^{YQg,YQi, A) of soft hadrons within sub-jet angle with respect to the energy flow in 
terms of the gluon jet multiplicity Nq [Yq , A) 



^l(>eo,>e,A) « —{C)Ao(.Yeo,Ye,X)N^iYe,X), 



(29) 



where A''^(l0, A) is given by the solution Eq. ([TT]) of the evolution equation (ID and {C)ao (^60 ^Y^e^X) 
(= {C)a(}{0) is average color current of partons forming the energy flux 



(C)A„(0=A^c(n)5o(0 + CF(n)^„(0 

+ 7o(i^e) \N,{nlnu)%{0 + CF ((nlnn)«^(0 + ri(n)^^(e) 



(30) 



70 (^0) 



N,[^{uln'u)%{0-ai{ulnu)%^{0 



G 



These color currents {C)ao are also presented in the appendix with numerical coefficients for = 3 
(see Eqs. (l40l ) and (|4TI)). The ratio of the gluon to the quark jet average multiplicity reads 



(31) 



For this quantity in the limit ^ 0o (C ^ 0), the appropriate ratio r = ^ = ro(l — ri7o — r27o) 
recovered. 

The color currents inside gluon and quark jets {C)c and (C)q are displayed in Fig.[2]for a particular 
jet virtuality £^©0 = 100 GeV in the LLA, MLLA and NMLLA. For simplicity we also use Qo = 
^QCD (A = 0) which is known to be a good approximation for hadron spectra ("limiting spectrum"). 
As expected from Eqs. (l30l ) and (l40l ). for ^ ©o — > 0), all curves for (C)G(^e) coincide at 
{C)c{YQg) = Nc, where Nc is the color factor corresponding to the splitting of one gluon (in this case 
the one initiating the jet) into two other gluons. In the same limit © —>■ ©0, the color current {C)q 
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Figure 2: Color currents {C)c (left) and {C)q (right) as a function of Yq = ln{EQ/Qo) for Yq < 
Yeo = InC-EBo/Qo) with £'9o = 100 GeV for different approximations (LLA, MLLA, NMLLA) with 
Qo = Aqcd = 230 MeV. 



in LLA approaches Cp, where Cp i& the color factor corresponding to the splitting of a quark into a 
quark and a gluon. The mismatch at = 0o between the MLLA and the LLA curves comes from the 
7o-expansion of r = Nq/Nq and equals Cpri'yo, while the difference between the LLA curve and the 
NMLLA curve is increased by the O(7o) correction and reaches Cpiri'^o + ^27o)- 
The role of 0(70) and O(7o) corrections to the LLA color current {C)ao in Eq. (l30l ) is quite important 
and the expansion in 70 is seen in Fig. [2] to oscillate if higher order terms are incorporated. Indeed, 
the terms 7Q(nln* u) oc (—70)* change sign as one goes from the LLA (i = 0) to MLLA (i = 1) and 
to the NMLLA (i = 2). Because of the running coupling the 70-expansion (70 ~ ^/al) converges at 
high energies but may diverge at low energies. At present energies (70 ~ 0.5 at LEP) the 70 expansion 
is converging rather slowly and one may ask for yet higher order terms. Such calculations would also 
require higher order terms in the expansion of the ratio r and its derivatives which are not available. 
As an exercise, we calculated 6 terms beyond LLA for the gluon current numerically while keeping the 
expansion of r at second order in 79. Then the solutions keep oscillating with a new solution falling in 
between the two previous ones. The trend is towards a solution closer to NMLLA than to MLLA. 

Starting with the isolated jet at Gq (^60) the sub-jets at smaller angles Q (Yq) may evolve from interme- 
diate partons of different color. Therefore, one expects that sub-jets in a gluon jet have lower multiplicity 
than isolated gluon jets at the same angle (Ye) and the opposite for a quark jet. This behavior is in fact 
bom out by the MLLA and NMLLA curves for {C)g which show the same trends as the LLA term: it is 
increasing with Yq at the higher energies. On the other hand, the quark color current (C)q is decreasing 
in LLA and in the NMLLA while for the intermediate approximation MLLA this behavior is not yet 
reached. Asymptotically (70 0) all curves have to approach either {C)g = Nc or (C)q = Cp- 

3.3 Results on multiplicities at small angles 

Next we study the consequences for the behavior of multiplicities at full angles Gq ~ 1 and for sub-jets 
at reduced angles G. Results for multiplicities using our formulae for color currents Eqs. (l29l ) and (l30l ) 
are shown in Figs. [3] and |4] for MLLA (left panels) and NMLLA (right panels). The full lines show the 
particle multiplicity A^^o for isolated Aq jet as function of y = Yq^ = ln(ii^Go/(5o)> i-S- as function 
of the jet energy E for the gluon jet in Fig. |3]and for the quark jet in Fig. |4l according to Eqs. © and 
(fTTI) . The hadronization constant has been taken as = 0.2 according to Also shown in these 
figures are the multiplicities for the reduced opening angle as function of Yq = hi{EQ / Qq) for 
different fixed energies E but variable sub-jet angles G. For G ^ Gq these curves approach the full 
multiplicities at y = leo- In the simple model with the scaling of multiplicities with EQ in Eq. ^ 
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Figure 3: Multiplicity in gluon jets as function of the variable Y = ln{E@ / Qq): the full line represents 
multiplicity of isolated jet Ng{EQq) at fixed opening angle G = 0o = 1 varying the energy E, the other 
curves represent sub-jet multiplicities Ng{E, Bq, ©) at fixed energies E = 50, 100, 1000 GeV varying 
the angles G; left panehin MLLA, right panel: in NMLLA {Qq = Aqcd = 0.23 GeV). 





Figure 4: Multiplicities as in Fig. [3l but for quark jets. 



all these curves would coincide; the existence of intermediate processes and non-tiivial color currents 
yields the predicted scale breaking. The mixing of quark and gluon jets will reduce the multiplicity in 
the gluon jet and increase it in the quark jet at sufficiently high energies. This property is reproduced in 
our calculations except for N^^^ in MLLA where the appropriate high energy regime is not yet reached. 

In case of the ratio r = Nq/Nq the effects from gluon and quark jets go in the same direction and yield 
the observable scale breaking effects of Fig. |5]for the small cone measurements. For the LEP energy 
range (Y ~ 5) the ratio for large jet opening angle ©o is r ~ 2 in MLLA and r ^ 1.8 in NMLLA. 
This is stiU larger than the experimentally observed r ~ 1.5 ITHH which indicates the slow convergence 
of the 70 ~ expansion of multiplicity. As emphasized in Sect. 2 the numerical treatment of the 
evolution equation removes this discrepancies largely. The effect from intermediate processes in the new 
definition of multiplicity amounts to about 20% in NMLLA at a reduced energy scale Y ^ Yq — 2. 



4 Conclusions 

We have studied a definition of jet axis which is based on a 2-particle correlation where the jet axis 
corresponds to the energy weighted direction of particles in a given cone ©q. In this way smaller angles 
(or transverse momenta) can be meaningfully determined. We studied especially the effects related to a 
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Figure 5: Multiplicity ratio Nq/Nq for jets and sub-jets as function of energy variable Y = \\i{EQ / Qq) 
as in Fig. [3] 



measurement of multiplicities where we have added MLLA and NMLLA corrections to the known LLA 
results and studied phenomenological consequences numerically. 

The main new effect from 2-particle correlations is the appearance of the "color current" in the result 
which reflects the possibility of intermediate quark-gluon processes. This effect can be seen by the scale 
breaking between jets and sub-jets at the same scale EQ but different energies E and opening angles 6. 
Typical effects are of the order of 20%. 

The expansion in ,Jal which is used in the calculation of color currents is rather slowly converging and 
MLLA in the considered energy range is not satisfactory whereas NMLLA shows the expected behav- 
ior at the qualitative level. For a more quantitative treatment higher order terms have to be calculated 
but this requires also a treatment of the multiplicity ratio r at higher orders which is not yet available. 
Alternatively, numerical methods and Monte Carlo calculations may be accessible. 

As an interesting application, one could learn from a single type of jet about properties of the other type 
by the study of sub-jets of lower energies (for example, from a high energy gluon jet at LHC about quark 
jets of lower energies). Other applications concern behavior of spectra at small transverse momenta or 
low energy gluon jets. 
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Appendix: Expressions for {u In* '^)ao^ the color currents {C)ao 

Setting = ,^(1) we obtain the MLLA correction from Eq. (l23l) for i = 1 and nj = 3 

(nlnu)g(^) = -0.875023 - 7.80247^ + e'^'^^^^'^^ (0.875023 - 1.65883^), (32) 

(nlntt)g(^) = 0.326411 - 4.38889^ + e-^-^^^^^^(-0.326411 + 1.65883^), (33) 

(nln'u)§(0 = 0.100287 - 7.80247^ + e-^-^^^^^^(-0.100287 + 2.94903^, (34) 

(uln'u)g(O = 0.875023 - 4.38889^ + e-^-^^^^^^(-0.875023 - 2.94903^; (35) 
for the NMLLA correction with i = 2 we also write the result forn j = 3 

(nln2n)g(^) = 7.64366e-^-^^^^^^(-1.0829 + ^(-0.146177 + ^ 

+ 95.1228(-0.0286673 + C)(0.443708 + 0, (36) 

(iiln2u)g(C) = 53.5066(-0.117852 +^)(0.159854 + ^) 

- 7.64366e-^-^^^^^^ (-0.744713 + 0(0.177083 + 0, (37) 

(nlii2n)§(0 = 95.1228(-0.0746644 + 0(0.239705 + 

- 13.5887e-^-^^^^^^ (-0.495132 + ^ (0.253032 + , (38) 

(nln^ n)g(0 = 1.20995 + 13.5887e-^-^^^^^« (-0.155024 + 0(0.574368 + 

+ 0-11-1292 + 53.50660- (39) 



Substituting the former results in (1301 ) and setting 

7o(i?e) = 7o(i^eo)e2^^'5o^ 7o(^0) = 7o'(i?0o)e^^^'^°« 
yields the color currents, evaluated again forn j = 3, 

(C)g(O = 2.4 + 0.6e-5-^5^5^« (40) 
+ 7o(£^eo) e~^-°^^^^^ (2.10105 + e^-^^^^''^(-2.10105 - 29.25930 - 2.76472^) 

+ 7o(-Eeo) (6.36971e^-^^^^^«(-1.39557 + O(-0.176893 + O 
+ 178.355 e^^(-0.0219197 + 0(0.402217 + 0) +0{j^), 



(C)q(O = 2.4- 1.06667 e-5-^5^5^« (41) 
+ 7o(Seo) (e~i °^55*^« (^-1.30969 + e5-5555^«(1.55636 - 29.25930 +4.91505^)) 
+ 7o (-EGo) (l78.355 ^(-0.0527303 + 0(0.183028 + - 11-3239 e3-44444C 



(-0.832991 + 0(0.247506 + 0) + C(7o)- 
After replacing by = 1) in CHI' one obtains the functions {C)c,Q(YeQ, Yq, A) involved in (|29l ) 
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